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Abstract 

In general or normal random matrix ensembles, the support of eigenvalues of large 
size matrices is a planar domain (or several domains) with a sharp boundary. This 
domain evolves under a change of parameters of the potential and of the size of 
matrices. The boundary of the support of eigenvalues is a real section of a complex 
curve. Algebro-geometrical properties of this curve encode physical properties of 
random matrix ensembles. This curve can be treated as a limit of a spectral curve 
which is canonically defined for models of finite matrices. We interpret the evolution 
of the eigenvalue distribution as a growth problem, and describe the growth in 
terms of evolution of the spectral curve. We discuss algebro-geometrical properties 
of the spectral curve and describe the wave functions (normalized characteristic 
polynomials) in terms of differentials on the curve. General formulae and emergence 
of the spectral curve are illustrated by three meaningful examples. 
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1 Introduction and Preliminciries 

Recently, random matrix theory has found new apphcations in growth prob- 
lems, where the evolution of the interface separating domains of different na- 
ture is a subject of interest. In some realizations, the growing domain is an 
aggregate of randomly deposited subunits. 

In general random matrix ensembles, complex eigenvalues usually occupy pla- 
nar domains in the complex plane. Their boundaries become sharp as the 
size of matrices, A^, goes to infinity. It appears that, for an important class of 
growth models, the aggregates evolve similarly to the support of eigenvalues of 
general matrix ensembles. In Refrs. [1,2,3], one of the most interesting classes 
of growth problems - Laplacian growth - has been linked to the evolution of 
normal random matrices. 

The interpretation of random matrix theory in terms of an aggregation pro- 
cess seems to be a productive approach in a number of different applications. 
Among them are the growth problems mentioned above, the semiclassical be- 
havior of electronic droplets in the Quantum Hall regime [2], and M = 1 
super symmetric Yang-Mills theory [4]. They reveal a relatively new, geomet- 
rical aspect of random matrix theory. From the mathematical point of view, 
the connection of random matrix ensembles with isomonodromic deformations 
of differential equations [5,6,7,8] seems to be especially interesting for these 
applications. 

It has been emphasized in recent papers [4] and [9] that the boundary of the 
domains of eigenvalues is to be thought of as a real section of a complex al- 
gebraic curve (a Riemann surface). This curve appears to be a fundamental 
object of random matrix theory. It encodes the most interesting physical prop- 
erties of matrix ensembles (see, e.g., [4,9,10,11] for the complex curves in one- 
and two- matrix models). 

In the Hermitian or unitary ensembles, the support of eigenvalues is a set of 
line intervals. Their boundaries are just points. In general matrix ensembles, 
eigenvalues are complex numbers, and their support is a planar domain. The 
boundaries are planar curves, evolving in a complicated and unstable man- 
ner. The notion of complex curve naturally links the planar geometry of the 
boundaries to the algebro-geometric properties of the matrix ensemble. 

In this paper we present general aspects of the normal matrix ensemble from 
this standpoint. One may canonically associate a complex curve to the matrix 
ensemble, not only in the large limit, but for any finite N as well. It is 
the spectral curve of the operator which generates recursion relations (the 
Lax operator), projected to a certain finite-dimensional subspace. In the large 
N limit (a semiclassical limit), this curve becomes the complex curve which 
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describes the support of eigenvalues. 

We introduce the spectral curve for normal matrix ensembles, and describe the 
evolution of the curve with respect to parameters of the statistical weight of 
the ensemble (deformation parameters) and the size of matrices (a parameter 
of growth) . In a forthcoming paper, we will prove that the semiclassical limit of 
the evolution (or deformation) equations is the universal Whitham hierarchy 
associated with the complex curve. The semiclassical wave function can be 
expressed through differentials on the curve. We also notice that the Whitham 
hierarchy is identical to the set of equations which describe Laplacian growth 
processes - unstable dynamics of an interface between two immiscible phases. 

Three meaningful examples illustrate the general formulae. 



1.1 Normal matrix ensemble 



A matrix M is called normal if it commutes with its Hermitian conjugate: 
[M, Mt] = 0, so that both M and can be diagonalizcd simultaneously. 
Eigenvalues of normal matrices are complex. The statistical weight of the 
normal matrix ensemble is given through a general potential W{M, M^) [12]: 

eb^mM,M^)d^^M). (1) 

Here /i is a parameter, and the measure of integration over normal matri- 
ces is induced by the flat metric on the space of all complex matrices dcM, 
where dcM — YiijdTZe M^jdZm Mij. Using the standard procedure, (see, e.g., 
[13]) one passes to the joint probability distribution of eigenvalues of normal 
matrices 2:1, . . . , ^jv, where N is size of the matrix: 

1 ^ 

-i-|A^(.)r (2) 

Here d^Zj = dxj dyj for Zj = Xj + iyj, Aj^{z) = det(^]~^)i<jj<7v = U.iij{zi-Zj) 



is the Vandermonde determinant, and 

1 r ^ 

r^ = ^/|A^(.)rn^^"^^^-'^-^^%- (3) 



is the normalization factor. This is the partition function of the matrix model 
(a T-function). 

A particularly important special case arises if the potential W has the form 

W ^-\z\^ + V{z) + V{z), (4) 
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Fig. 1. A support of eigenvalues consisting of four disconnected components. 




Fig. 2. The distribution of eigenvalues for the Gaussian potential. The droplet is an 
ellipse with quadrupole moment 2\t2\ and area irhN. 

where V{z) is a holomorphic function in a domain which includes the support 
of eigenvalues (see also a comment in the end of Sec. 1.6 about a proper 
definition of the ensemble with this potential). In this normal matrix 

ensemble gives the same distribution as a general complex matrix ensemble. 
A general complex matrix can be decomposed as M = U{Z + R)U'^ , where 
U and Z are unitary and diagonal matrices, respectively, and R is an upper 
triangular matrix. The distribution (2) holds for the elements of the diagonal 
matrix Z which are eigenvalues of M (Ref. [13]). Here we mostly focus on the 
special potential (4), and also assume that the field 



(a "vector potential", see Sec. 1.5) is a globally defined meromorphic function. 

1.2 Droplets of eigenvalues 

In a proper large limit {h 0, Nh fixed), the eigenvalues of matrices from 
the ensemble densely occupy a connected domain D in the complex plane, 
or, in general, several disconnected domains. This set (called the support of 
eigenvalues) has sharp edges (Fig. 1). We refer to the connected components 
Da of the domain D as droplets. 

In the case of algebraic domains (the definition follows) the eigenvalues are 
distributed with the density p = —-^AW, where A = 4:8^8^ is the 2-D Laplace 
operator [14]. For the potential (4) the density is uniform. The shape of the 
support of eigenvalues is a more involved subject. We discuss it below. For 
example, if the potential is Gaussian [15], 




(5) 



A{z) = 2t2Z, 



(6) 
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Fig. 3. The distribution of eigenvalues for the potential V{z) = —alog{l — z/P)—'yz. 
If the area of the droplet is ttHN, the complement of the droplet is an algebraic 
domain which can be conformally mapped to the exterior of the unit disk by the 
Joukowsky map (61). In this case, the cut of the Schwarz function shown in the 
figure (corresponding to a virtual droplet located on the unphysical sheet) shrinks 
to a double point. 




Fig. 4. The distribution of eigenvalues for the potential V{z) = tsz^. Cuts of the 
Schwarz function are shown. The cuts correspond to virtual droplets located on 
unphysical sheets. If the area of the droplet is nhN, the cuts shrink to double 
points. The boundary contour is a hypotrochoid. 

the domain is an ellipse (see Fig. 2). If ^4 has one simple pole, 

^«^-^-7 (T) 

the droplet (under certain conditions discussed below) has the profile of an 
aircraft wing given by the Joukowsky map (Fig. 3). li A has one double pole 
(say, at infinity), 

A{z) = 3t3z\ (8) 

the droplet is a hypotrochoid (Fig. 4). These are the three examples that we 
will discuss below in detail. If A has two or more simple poles, there may be 
more than one droplet. 

If A{z) is a more complicated function, the domain of eigenvalues develops an 
unstable fingering pattern, similar to the one in Fig. 5. 

Boundary components of the droplets form a real section of a complex curve. 
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Fig. 5. A grown fingering pattern observed in a radial Hele-Shaw cell. Air is inserted 
under pressure to the cell filled by silicon oil [16]. 

1.3 Semiclassical complex curve 

Jumping ahead, we describe the construction of the semiclassical complex 
curve which emerges in the normal matrix ensemble. Let us represent the 
boundary of the domain as a real curve F{x, y) = 0. If the vector potential 
A{z) is a meromorphic function (we always assume that this is the case), the 
function F can be chosen to be an irreducible polynomial. Then we rewrite it 
in holomorphic coordinates as 



and treat z and z as independent complex coordinates z, z. The equation 
f{z,z) = defines a complex curve. This curve is a finite-sheet covering of 
the z-plane. The single-valued function z{z) on the curve is a multivalued 
function on the 2;-plane. Making cuts, one can fix single-valued branches of 
this function. 

The boundary of the domain is a section of the curve by the plane where z is 
complex conjugate of z. It belongs to a particular sheet (we call it the physical 
sheet) of the covering. It appears that physical properties of the ensemble are 
determined not only by the physical sheet, but rather by the entire algebraic 
covering, including all the sheets other than physical. The Riemann surface 
for Example (7) is presented in Fig. 10. 

1.4 Growth and deformation parameters. Evolution of the curve 

Similarly to the Hermitian matrix ensemble [4,9,10,11], the complex curve of 
the normal matrix ensemble is characterized by the potential W, or by the 
"vector potential" A{z), and by a set of gf-|- 1 integers (not necessarily pos- 
itive), where g is genus of the curve. The integers are subject to the constraint 




(9) 
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Yl^=o — N ■ We will discuss the meaning of these numbers later on. For now, 
we only mention that if they are all positive, then they are proportional to 
the areas of the droplets of uniformly distributed eigenvalues. Every droplet 
contains Ua eigenvalues, so a quantum of area is uTi per particle. 

As one varies the potential and the filling factors Va-, the curve and the inter- 
face bounding the droplets evolve. Parameters of the potential (for example, 
poles and residues of the meromorphic function (5)) and filling numbers are 
deformation parameters and parameters of growth. They are coordinates in 
the moduli space of the complex curves. 

An infinitesimal variation of the potential generates correlation functions of 
the ensemble (Ref. [14]). For example, irreducible correlation functions of re- 
solvents (holomorphic currents) J{z) — tr = Y.i are generated by 
variations of A(z). In Ref. [14] we showed that the correlation functions are 
expressed through algebro-geometrical properties of the boundary of the do- 
main. 

During the evolution, the genus of the complex curve may change. This cor- 
responds to a coalescence of droplets or a droplet breakup, or even more com- 
plicated degenerations, when a droplet shrinks to a point and simultaneously 
merges with another droplet. Close to degeneracy points (also called critical 
points), matrix ensembles have universal scaling behavior. We will present a 
study of this problem in a future paper. 

A particularly interesting process is growth. In this case, one changes the total 
number of eigenvalues N and keeps the potential fixed. Below we refer to the 
normalized total number of particles 

t = hN 

as time of growth. In the algebraic case, the time t is the normalized (modulo 
tt) area of the droplets. While increases, the new eigenvalues aggregate at 
the boundary of the existing droplets, so that the domain of eigenvalues (an 
aggregate) grows (Fig. 6). Growth of a hypotrochoid is depicted in Fig. 7. 

There are several important physical problems where this kind of evolution 
occurs. We list them in historical order. 

The first one is the celebrated Hele-Shaw problem - the most studied example 
of Laplacian growth. The Hele-Shaw problem describes the non-equilibrium 
dynamics of a planar interface between two immiscible fluids conflned to a 
thin 2-D cell. The second one is the matrix model description of 2-D quantum 
gravity. The third one is the evolution of energy levels in mesoscopic systems. 
A similar problem also emerges in isomonodromic deformation of differential 
equations. 
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Fig. 6. The area of a droplet grows with the rate irfi per eigenvalue. A new eigenvalue 
aggregates at the boundary of the droplet. The shaded area depicts the support of 
the semiclassical wave function. 




Fig. 7. Hypotrochoid grows until it reaches a critical point. 

The last three problems and their connection with the matrix ensemble are 
well reflected in the literature (see e.g., [7,8,17] for reviews). The connection 
of matrix ensembles to the Hele-Shaw problem is relatively new. It nicely 
illustrates the growth and gives a transparent hydrodynamic interpretation 
to abstract objects of algebraic geometry of curves and complex analysis on 
Riemann surfaces. We find that a brief description of the Hele-Shaw problem 
in this paper is in order (see Sec. 7). 



1.5 Normal matrix ensemble and quantum Hall effect 

A useful interpretation of the Coulomb gas distribution (2) is a coherent state 
of relativistic electrons in the Quantum Hall regime [2]. In this case, the elec- 
trons are situated in the plane in a strong, not necessarily uniform, magnetic 
field B{z) — —^AW, and fully occupy the lowest energy level. The exact 
A^-particle wave function, defined up to a phase, is 

^n{zi, ...,z^)^ -^^Ar,{z) e^U ^^(--^.■). (lo) 
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The joint probability distribution (2) is then equal to |^(^i, . . . , zn)\'^- It is the 
probability to find electrons at the points Zi. In the semiclassical regime, the 
wave function (10) describes an incompressible electronic droplet with sharp 
edges. The function A{z) (5) can be thought of as a holomorphic component of 
the vector potential generated by magnetic impurities located far away from 
the electronic droplet. 

In this language, the growth problem translates into the evolution of a semi- 
classical electronic droplet under a change of magnetic field or chemical po- 
tential. 



1.6 Orthogonal polynomials as a measure of growth 



Let the number of eigenvalues (particles) increase while the potential stays 
fixed. If the support of eigenvalues is simply-connected, its area grows as HN. 
One can describe the evolution of the domain through the density of particles 

p^{z)^N j \-^N{z,Zi,Z2,...,ZN~l)\'^(fzi...(fzN-l, (H) 

where ^^at is given by (10). 

We introduce a set of orthonormal one-particle functions on the complex plane 
as matrix elements of transitions between N and [N + l)-particle states: 



i'N{z] 



= = j *JV+l(^, Zi,Z2,..., zn)'^n{zu -22, ... , ZN)d^Zi . . . Zn (12) 



Then the rate of the density change is 

pN+i{z) - pn{z) = \^n{z)\'^- (13) 

The proof of this formula is based on the representation of the tpn through 
holomorphic biorthogonal polynomials Pn{z). Up to a phase 



M^) = e^^(^'^")P„(^), Pn{z) = A^z^ + . . . (14) 

The polynomials Pn(^) are biorthogonal on the complex plane with the weight 

y e'^l''P^{z)P^{z)d\ = (15) 

The proof of these formulae is standard in the theory of orthogonal polyno- 
mials. Extension to the biorthogonal case adds no difficulties. 

A physical interpretation of the wave function is clear in the QHE setup. 
The eigenvalues Zi are the positions of electrons in a strong magnetic field. 
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Then |'?/'Ar(^)p is a probability of adding an additional electron to an aggregate 
of electrons at the point z. Since the droplet is incompressible, it is only 
possible to add the particle to the boundary of the droplet. We will see that 
the wave function in a semiclassical limit is indeed localized on the boundary. 

We note that, with the choice of potential (4), the integral representation (15) 
has only a formal meaning, since the integral diverges unless the potential is 
Gaussian. A proper definition of the wave functions goes through recursive 
relations (17, 18) which follow from the integral representation. The same 
comment applies to the r-function (3). The wave function is not normalized 
everywhere in the complex plane. It may diverge at the poles of the vector 
potential field. 

Below we describe the evolution and semiclassical behavior of the wave func- 
tion in the region of the complex plane close to the boundary of physical 
droplets. 



2 Equations for the wave functions and the spectral curve 

In this section we specify the potential to be of the form (4). It is convenient 
to modify the exponential factor of the wave function. Namely, we define 

V'n(^) = e-^+^^Wp„(z), and Xn(^) = ei^Wp„(z), (16) 

where the holomorphic functions Xn(-2) are orthonormal in the complex plane 
with the weight e~'^' Like traditional orthogonal polynomials, the biorthog- 
onal polynomials P„ (and the corresponding wave functions) obey a set of 
differential equations with respect to the argument and recurrence relations 
with respect to the degree n. Similar equations for two-matrix models are 
discussed in numerous papers (see, e.g., [18]). 

We introduce the L-operator (the Lax operator) as multiplication by z in the 
basis Xn- 

LnmXm{z) = ZXn{,z) (17) 

(summation over repeated indices is implied). Obviously, L is a lower triangu- 
lar matrix with one adjacent upper diagonal, Lnm = Oasm>n-|-l. Similarly, 
the differentiation dz is represented by an upper triangular matrix with one 
adjacent lower diagonal. Integrating by parts the matrix elements of the 9^, 
one finds: 

{L'^)nniXm = fld^Xn, (18) 

where is the Hermitian conjugate operator. 

The matrix elements of are {L^)nm — -^mn — ^{Lnm)+ 1 Pm(z)dzPn{z)(Pz, 
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where the last term is a lower triangular matrix. The latter can be writ- 
ten through negative powers of the Lax operator. Writing dz log Pn{z) — 
7 + Z]a;>i one represents L"^ in the form 

^A{L) + {nn)L-^ + ^v^''^L-\ (19) 

k>l 

where f*^^-* and (hn) are diagonal matrices with elements vlj'^ and (hn). The 
coefficients v^^ are determined by the condition that lower triangular matrix 
elements of A{Lnm) are cancelled. 

In order to emphasize the structure of the operator L, we write it in the basis 
of the shift operator ^ w such that wfn = fn+i^ for any sequence /„. Acting 
on the wave function, we have: 

WXn = Xn+l- 

In the n-representation, the operators L, acquire the form 

L = rnW + Y, ul^^w-\ = w-'rn + ^ w''u^J^\ (20) 

fe>0 fe>0 

Clearly, acting on x„, we have the commutation relation ("the string equa- 
tion") 

[L, Lt] = h. (21) 
This is the compatibility condition of Eqs. (17) and (18). 

Equations (20) and (21) completely determine the coefficients and u!^\ 

The first one connects the coefficients to the parameters of the potential. The 
second equation is used to determine how the coefficients r„ and u!;!^'^ 
evolve with n. In particular, the diagonal part of it reads 

nh - rl-j: j:\uUf . (22) 
fc>ip=i 

Moreover, we note that all the coefficients can be expressed through the r- 
function (3) and its derivatives with respect to parameters of the potential. 
This representation is particularly simple for r„: = 



^ The shift operator w has no inverse. Below iv^^ is understood as a shift to the 
left defined as w''^w = 1. Same is applied to the operator L^^. To avoid a possible 
confusion, we emphasize that although Xn is a right-hand eigenvector of L, it is not 
a right-hand eigenvector of L"^. 
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2.1 Finite dimensional reductions 



If the vector potential A{z) is a rational function, the coefficients it^'^) are not 
all independent. The number of independent coefficients equals the number 
of independent parameters of the potential. For example, if the holomorphic 
part of the potential, V{z), is a polynomial of degree d, the series (20) are 
truncated at k = d — 1. 

In this case the semi-infinite system of linear equations (18) and the recurrence 
relations (17) can be cast in the form of a set of finite dimensional equations 
whose coefficients are rational functions of z, one system for every n > 0. The 
system of differential equations generahzes the Cristoffel-Daurboux second 
order differential equation valid for orthogonal polynomials. This fact has been 
observed in recent papers [19,20] for biorthogonal polynomials emerging in the 
Hermitian two-matrix model with a polynomial potential. It is applicable to 
our case (holomorphic biorthogonal polynomials) as well. 

In a more general case, when A{z) is a general rational function with d — 1 
poles (counting multiplicities), the series (20) is not truncated. However, L 
can be represented as a "ratio" , 

L = K^^K2 = MaMf \ (23) 

where the operators i^i,2, Afi,2 are polynomials in w: 

d-2 d-l 

Ki = w''-' + J2 ^n^^' > ^2 = rr^iw'' + ^ B^^^w^ (24) 

j=0 j=0 
d-2 d-l 

3=0 3=0 

These operators obey the relation 

K1M2 = K2M1. (26) 

It can be proven that the pair of operators Mi^2 is uniquely determined by 
Ki^2 and vice versa. We note that the reduction (23) is a difference analog of 
the "rational" reductions of the Kadomtsev-Petviashvili integrable hierarchy 
considered in [21]. 

The linear problems (17), (18) acquire the form 

{K2X)n = Z {Kix)n , (Mlx)n = hd,{Mlx)n. (27) 

These equations are of finite order (namely, of order d), i.e., they connect 
values of x„ on d -|- 1 subsequent sites of the lattice. 
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The semi- infinite set {xo, Xi? ■ ■ ■ } is then a "bundle" of d-dimensional vectors 

Xin) = iXn, Xn+l, • • • , Xn+d-lY 

(the index t means transposition, so x is a column vector). The dimension of 
the vector is the number of poles of A{z) plus one. Each vector obeys a closed 
(i-dimensional linear differential equation 

hd,x{n)^Cn{z)x{n), (28) 

where the d x d matrix Cn is a "projection" of the operator onto the n-th 
d-dimensional space. Matrix elements of the Cn are rational functions of z 
having the same poles as A{z) and also a pole at the point A{oo). (If A{z) is 
a polynomial, all these poles accumulate to a multiple pole at infinity). 

We briefly describe the procedure of constructing the finite dimensional matrix 
differential equation. We use the first linear problem in (27) to represent the 
shift operator as a d x d matrix W„(^) with ^-dependent coefficients: 

Wn{z)x{n)^x{n+1). (29) 

This is nothing else than rewriting the scalar linear problem in the matrix 
form. Then the matrix Wn{z) is to be substituted into the second equation of 
(27) to determine Cn{z) (examples follow). The entries of Wn{z) and Cn{z) 
obey the Schlesinger equation, which follows from compatibility of (28) and 
(29): 

hd.Wn = Cn+l'^n ' ^^nC^n- (30) 

This procedure has been realized explicitly for polynomial potentials in re- 
cent papers [19,20]. We will work it out in detail for our three examples: 
x(n) = (x„, Xn+i)* for the ellipse (6) and the aircraft wing (7) and x(n) = 
(Xn, Xn+i,Xn+2)* for the hypotrochoid (8). 

2.2 Spectral curve 

According to the general theory of hnear differential equations, the semiclassi- 
cal (WKB) asymptotics of solutions to Eq. (28), as ^ — > 0, is found by solving 
the eigenvalue problem for the matrix Cn{z) (see, e.g., [22]) . More precisely, 
the basic object of the WKB approach is the spectral curve [23] of the ma- 
trix Cn, which is defined, for every integer n > 0, by the secular equation 
det{jCniz) — z) — (here z means z ■ 1, where 1 is the unit d x d matrix). It 
is clear that the left hand side of the secular equation is a polynomial in z of 
degree d. We define the spectral curve by an equivalent equation 

fn{z,z) = a{z)det{Cn{z) - z) = 0, (31) 
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where the factor a{z) is added to make fn{z, z) a polynomial in z as well. The 
factor a{z) then has zeros at the points where poles of the matrix function 
C{z) are located. It does not depend on n. We will soon see that the degree 
of the polynomial a{z) is equal to d. Assume that all poles of A{z) are simple, 
then zeros of the a{z) are just the d—1 poles of A{z) and another simple zero 
at the point A{oo). Therefore, we conclude that the matrix Cn{z) is rather 
special. For a general d x d matrix function with the same d poles, the factor 
a{z) would be of degree d^. 

Note that the matrix jCni^) ~ ^ enters the differential equation 

hd,\l{n)\^ = l{n){£n{z) - z)l{n) (32) 

^|2 

for the squared amplitude = ijj\n)ijj{n) = e «~|x(n)p of the vectors 

il^iji) built from the orthonormal wave functions (14). 

The equation of the curve can be interpreted as a "resultant" of the non- 
commutative polynomials K2 — zKi and — zM\ (cf. [19]). Indeed, the 
point (2;, z) belongs to the curve if and only if the linear system 

{K2c)k = z{i^\c)k n — d <k <n — \ 

(33) 

{Mlc)k = z{M[c)k n<k<n + d-l 

has non-trivial solutions. The system contains 2d equations for 2d variables 
Cn-d , ■ ■ ■ , Cn+d-i- Vanishing of the 2d x 2d determinant yields the equation of 
the spectral curve. Below we use this method to find the equation of the curve 
in the examples. It appears to be much easier than the determination of the 
matrix Cn{z)- 

The spectral curve (31) possesses an important property: it admits an antiholo- 
morphic involution. In the coordinates z, z the involution reads {z, z) 1-^ {z, z). 
This simply means that the secular equation det{Cn{z) — z) — for the matrix 
£„(5) = Cn{z) defines the same curve. Therefore, the polynomial /„ takes real 
values for 5 = f : 

U{z,z) = U^. (34) 
Points of the real section of the curve {z = z) are fixed points of the involution. 

The curve (31) was discussed in recent papers [19,20] in the context of Her- 
mitian two-matrix models with polynomial potentials. The dual realizations 
of the curve pointed out in [19] correspond to the antiholomorphic involution 
in our case. The involution can be proven along the lines of these works. The 
proof is rather technical and we omit it, restricting ourselves to the examples 
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below. We simply note that the involution relies on the fact that the squared 
modulus of the wave function is real. 



2.3 Schwarz function 



The polynomial fn{z, z) can be factorized in two ways: 




(35) 



Uzrz) = a{z){z-S^^\z))...{z-S^^\z)l 



(36) 



where S^\z) are eigenvalues of the matrix Cn{z). One may understand them 
as different branches of a multivalued function S{z) (respectively, S{z)) on 
the plane (here we do not indicate the dependence on n, for simplicity of the 
notation). It then follows that S{z) and S{z) are mutually inverse functions: 



An algebraic function with this property is called the Schwarz function. By the 
equation f{z,S{z)) = 0, it defines a complex curve with an antiholomorphic 
involution. An upper bound for genus of this curve is g = {d — 1)^, where d 
is the number of branches of the Schwarz function. The real section of this 
curve is a set of all fixed points of the involution. It consists of a number of 
contours on the plane (and possibly a number of isolated points, if the curve is 
not smooth). The structure of this set is known to be complicated. Depending 
on coefficients of the polynomial, the number of disconnected contours in the 
real section may vary from to (? + 1. If the contours divide the complex curve 
into two disconnected "halves", or sides (related by the involution), then the 
curve can be realized as the Schottky double [24,25,26] of one of these sides. 
Each side is a Riemann surface with a boundary. 

We will discuss general properties of the Schwarz function and the Schottky 
double in Sec. 4. 

Let us come back to equation (28). It has d independent solutions. They are 
functions on the spectral curve. One of them is a physical solution correspond- 
ing to biorthogonal polynomials. The physical solution defines the "physical 
sheet" of the curve. 

The Schwarz function on the physical sheet is a particular root, say S^\z), of 
the polynomial fn{z, z) (see (35)). It follows from (19) that this root is selected 
by the requirement that it has the same poles and residues as the potential A. 



S{S{z)) 



— z. 



(37) 
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A formal ^ semiclassical asymptote of the solution, in the leading order in h, 
is 

Xn-ey^'-\ (38) 

Here 

The differential dQ'^^^ is a physical branch of the generating differential on the 
curve (see below). 

The semiclassical asymptotics is discussed in more details in Sec. 6. Here wc 
notice that the amplitude of the wave function ipn{z) ( not Xn(-2)!) peaks at 
the solution of the equation 

S^'\z)^z, (39) 

where now z = x + iy and z = x — iy arc complex conjugated coordinates in 
the plane. Solutions to this equation describe those contours of a real section 
of the complex curve which belong to the physical sheet. It is a set of closed 
planar curves, as shown in Fig. 1. These curves are boundaries of semiclassical 
droplets. Evolution and growth of the droplets translate into evolution of the 
complex curve. 

2.4 Deformation parameters 

The essential information about the spectral curve (which, in what follows, 
we will call quantum curve, unless the semiclassical limit is considered) is 
contained in two conditions: 

• the antiholomorphic involution, and 

• poles and residues of the Schwarz function on one selected sheet (the physical 
sheet) are given by the vector potential field A{z). 

These requirements determine all but g coefficients of the polynomial fn{z, z), 
where g is genus of the curve. The remaining g coefficients do not depend on 
the potential (deformation parameters). 

The Bohr-Sommerfield quantization condition for the semiclassical wave func- 
tion (38) requires that the integral over all the cuts (or boundaries of the 
droplets) are integer multiples of 27ri. Also, the integral around the point at 
infinity, counting the degree of the polynomial, must be 2mN. The integrals 
over the cuts transform into the integrals over a-cycles of the curve. Later we 

^ This formal expression ignores the Stokes phenomenon. It is valid only around 
boundaries of physical droplets. 
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will interpret the a-cycles as boundaries of the droplets. Therefore, we have 
the set of integers 



The d-dimensional linear differential equations (28) give an interpretation of 
the potential data as deformation parameters. Varying the poles and residues 
of A(z), one does not change the monodromy of solutions, given by A'^ and i^a- 
They are similar to deformation parameters of the theory of isomonodromic 
deformations of the system of linear differential equations [27] (see also Ref. 
[8] for applications of the theory of isomonodromic deformations to Hermitian 
matrix models and orthogonal polynomials in one variable, and references 
therein). A connection between isomonodromic deformations and Whitham 
equations, which seems to be relevant to our discussion, was pointed out in 



In the spirit of the theory of isomonodromic deformations of linear differential 
equations, we conjecture that a specific value of N, a potential V{z), and a 
set of g numbers (40) form a complete data set which uniquely determines the 
coefficients u!^^ of the recurrence equations, and the solution of the differ- 
ential equation (28). We do not know, however, whether the requirement that 
the physical solutions are polynomials gives a restriction on possible values 
of the integers It is likely that a general set of Ua (positive and negative) 
corresponds to more general solutions of (28) rather than polynomials. 

We also conjecture that, for physical values of the parameters of the potential 
in general position, the quantum curve is smooth (i.e., non-degenerate), and 
has the maximal possible number of the disconnected closed contours ("real 
ovals") in the real section. In the latter case, the involution is known to di- 
vide the complex curve into two "halves" , thus making the Schottky double 
construction applicable. 

Before discussing general properties of the curve, we illustrate its appearance 
in the semi classical limit. 



3 Semiclassical curve 

By semiclassical case we mean the limit — oo, ^ — 0, while the time t = 
hN, and the potential remain fixed. In this case, all eigenvalues arc distributed 
continuously within droplets with sharp boundaries. The boundary of the 
droplets is a real section of the classical curve. 




a 



(40) 



Ref. [28]. 
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3.1 Saddle point equation and Riemann-Hilbert problem 

In the semiclassical limit, the distribution (2) peaks at the minimum of the 
function W{z) — (fi{z), where 

N 

cp[z) = -^^logl^ - Zi\'^ 

i=l 

is the Newton potential, logarithmic in 2-D, of the domain (the support of 
eigenvalues). In the semiclassical limit, the potential (p is continuous. One 
writes the extremum condition in the form 

d:,(ip-W) ^d^(ip-W) ^0, at Zi. (41) 

The points where this equations holds are the most probable positions of 
eigenvalues Zi. 

In addition, if one assumes that Zi have a continuous distribution (a strong 
assumption), i.e., that the density 

N -y 

Pm{z) ^hY,S{z-Zi)^- — (42) 
i=i 47r 

is a smooth positive function in some domain D, and zero outside this domain, 
then it follows from (41) that the density is —^AW inside the domain, if 
—AW is not singular and positive inside the support of eigenvalues. 

The shape of the domain is found in two steps. First, one finds the potential 
(p. Then, having the potential, one has to restore the shape of the domain. 
The second step is the inverse potential problem. 

Under this assumption, (/? is a harmonic function in the exterior of the support 
of eigenvalues. Then equation (41) reads as a Riemann-Hilbert problem: 

(i) The boundary value of the analytic function dz(p in the exterior of the 
domain D is dzW; 

(ii) The function d(p behaves at infinity as —hNz~^ — —t/z (we assume that 
there are no eigenvalues around infinity). 

These conditions uniquely specify the domain if it is simply connected. In this 
case, irhN is the area of the domain. Multiply-connected domains are also 
possible solutions. The number of disconnected droplets cannot exceed the 
number of poles of A{z) plus one. To find them one needs additional data. 
For example, one may fix the filling factors - the number of eigenvalues that 
each droplet contains. Obviously, all the filling factors are positive and their 
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sum is N. Domains constructed this way, provided A{z) is a rational function, 
belong to a special class of algebraic domains [29] . 

An important property of the Riemann-Hilbert problem is that the function 
d(p{z), with z inside the domain, does not depend on N, but rather on the 
potential. The former does not change during the growth. 

At this point, we specify the potential to be of the form (4). Then the boundary 
value of the analytic function d(p{z) is dW — —z + A{z). The density of 
eigenvalues is uniform, and the filling factors fix the areas of the droplets to 
be Trhua- 

The potential ip{z) can be constructed with the help of the Schwarz function. 

3.2 The Schwarz function and the inverse potential problem 

Given a closed (in general a multiply-connected) domain D in the plane, the 
Schwarz function S{z) is an analytic function in some neighborhood of the 
boundary of the domain with the value z on the boundary [24] : 

S{z) — z on the contour. (43) 

The involution property of the Schwarz function follows from the definition 
(cf. (37)): 

SiSiz)) = z. 

Let us represent the Schwarz function in the form S{z) = S-^{z) + S-{z), 
where S+ (respectively, -S"-) is analytic inside (respectively, outside) D, with 
the condition that 5'_(oo) = 0. Choose S+{z) to be equal to A{z), 

A{z)^^i zeD. (44) 

^ ' 2mJdDZ-C 

Then the function S^{z) = S{z) — A{z), being analytic outside D, has the 
boundary value z — A{z). We conclude that A{z) — S{z) = dip{z). 

Therefore, the problem of the previous section reads: given a potential A{z), 
find the Schwarz function 

S{z) = A{z) + S4z) (45) 

whose poles and residues outside of the support of eigenvalues are given by 
the poles and residues of A{z) and such that S-{z) — > hN/z as z ^ oo. If 
the Schwarz function is known, Eq. (43) determines the domain. Inside the 
droplet S-{z), and therefore, S{z) have cuts (Fig. 8). 
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Fig. 8. The boundary of a droplet and a cut of the Schwarz function located inside 
the droplet. 

A note is in order. In Sec. 3.1 we assumed that the density of eigenvalues is 
a continuous function everywhere except for the boundary of As a conse- 
quence, is a harmonic function, dip and S- [z) are analytic functions outside 
D, and the Schwarz function outside the droplets has the same poles as A{z). 
As a result, the complement to the support of eigenvalues appears to be an 
algebraic domain. 

It seems that this assumption is too restrictive. Moreover, the definition of 
the curve through the linear differential equations for the wave function of 
Sec. 2.2 does not forbid the physical branch of the Schwarz function to have 
branching points and cuts outside physical droplets, as shown in Figs. 3, 4. 
The semiclassical wave function still peaks on the real section of the curve 
(39), regardless of whether there are cuts outside the droplets, or if those cuts 
shrink to points. In other words, a solution with cuts outside physical droplets 
still has a local maximum on the boundary of the physical droplets and it 
is stable. It becomes unstable, however, in a vicinity of cuts along a normal 
direction to a cut. 

These arguments suggest that one can safely allow dip{z) to have branch- 
ing points, i.e., the Schwarz function is allowed to have cuts outside physical 
droplets. Then we seek for solutions of the Riemann-Hilbert problem of Sec. 
3.1 in a wider class of functions. This leads to a wider class of domains (not 
algebraic). However, in this case, one is no longer able to treat eigenvalues as 
continuously distributed. We will come back to this issue in the next section. 



4 The Schwarz function and its Riemann surface. The Schottky 
double. 

The Schwarz function describes more than just the boundary of clusters of 
eigenvalues. Together with other sheets it defines a Riemann surface. If the 
potential A[z) is meromorphic, the Schwarz function is an algebraic function. 
It satisfies a polynomial equation f{z, S{z)) = 0. 



20 




Fig. 9. The Schottky double. A Riemann surface with boundaries along the droplets 
(a front side) is glued to its mirror image (a back side). 

The function f{z,z), where z and z are treated as two independent complex 
arguments, defines a Riemann surface with antiholomorphic involution (34). 
If the involution divides the surface into two disconnected parts, as explained 
above, the Riemann surface is the Schottky double [25,26] of one of these parts. 

There are two complementary ways to describe this surface. One is through 
the algebraic covering (35, 36). Among d sheets we distinguish a physical sheet. 
The physical sheet is selected by the condition that the differential S{z)dz has 
the same poles and residues as the differential of the potential A{z)dz (as in 
(45)). It may happen that the condition z — S^'^^z) defines a planar curve (or 
several curves, or a set of isolated points) for branches other than the physical 
one. Wc refer to the interior of these planar curves as virtual (or unphysical) 
droplets situated on sheets other than physical. 

Another way emphasizes the antiholomorphic involution. Consider a mero- 
morphic function h{z) defined on a Riemann surface with boundaries. We call 
this surface the front side. The Schwarz refiection principle extends any mero- 
morphic function on the front side to a mcromorphic function on the Riemann 
surface without boundaries. This is done by adding another copy of the Rie- 
mann surface with boundaries (a back side), glued to the front side along the 
boundaries. Fig 9. The value of the function h on the mirror point on the back 
side is h{S{z)). The copies are glued along the boundaries: h{z) = h{S{z)) 
if the point z belongs to the boundary. The same extension rule applies to 
differentials. Having a meromorphic differential h{z)dz on the front side, one 
extends it to a meromorphic differential h{S{z))dS{z) on the back side. 

This definition can be applied to the Schwarz function itself. We say that the 
Schwarz function on the double is S{z) if the point is on the front side, and z 
if the point belongs to the back side (here we understand S{z) as a function 
defined on the complex curve, not just on the physical sheet). 
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Fig. 10. Physical and unphysical droplets on a torus. The physical sheet (shaded) 
meets the unphysical sheet along the cuts. The cut situated inside the unphysical 
droplet appears on the physical sheet. The boundaries of the droplets (physical and 
virtual) belong to different sheets. This torus is the Riemann surface corresponding 
to the ensemble with the potential V{z) = — alog(l — z/ (3) — ^z. 

The number of sheets of the curve is the number of poles (counted with their 
multiplicity) of the function A{z) plus one. Indeed, according to (45) poles of 
A are poles of the Schwarz function on the front side of the double. On the 
back side, there is also a pole at infinity. Since S{z = oo) = A{oo), we have 
S{z = A{oo)) = oo. Therefore, the factor a{z) is a polynomial with zeros at 
the poles of A{z) and at A{oo), and 

d = number of sheets = number of poles of A + 1. 

The front and back sides meet at planar curves z = S{z). These curves are 
boundaries of the droplets. We repeat that not all droplets are physical. Some 
of them may belong to unphysical sheets. Fig. 10. 

Boundaries of droplets, physical and virtual, form a subset of the a-cycles on 
the curve. Their number cannot exceed the genus of the curve plus one: 

number of droplets < g + I. 

The sheets meet along cuts located inside droplets. The cuts that belong to 
physical droplets show up on unphysical sheets. On the other hand, some cuts 
show up on the physical sheet (Fig 10). They correspond to droplets situated 
on unphysical sheets. 

The Riemann-Hurwitz theorem computes the genus of the curve as 
g = half the number of branching points — d + 1. 

With the help of the Stokes formula, the numbers {i^a} are identified with areas 
of the droplets: |i^a| = 2th -foa -^^^ ^ nondegenerate curve, these numbers 
are not necessarily positive. Negative numbers correspond to droplets located 
on unphysical sheets. In this case, do not correspond to the number of 
eigenvalues located inside each droplet, as it is the case for algebraic domains, 
when all filling numbers are positive. 
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Fig. 11. Degenerate torus corresponds to the algebraic domain for the Joukowsky 
map. 

A comment is in order. In the Hermitian or unitary one-matrix ensembles, the 
construction of the Riemann surface is different (Refrs. [4,10,11,30]). In the 
Hermitian random matrix ensemble, for example, eigenvalues are real. The 
support of eigenvalues is a set of disjoint intervals on the real line. In this 
case, branch cuts are identical to the support of eigenvalues. The Riemann 
surface is constructed by gluing g sheets along g + 1 cuts [30] . Conversely, in 
the Schottky double the cuts generally do not touch boundaries of the droplets 
(Fig. 8). Complex curves for the Hermitian or unitary ensembles are simpler. 
They are always hyperelliptic. 



Degeneration of the spectral curve 

Degeneration of the complex curve gives the most interesting physical aspects 
of growth. There are several levels of degeneration. We briefly discuss them 
below. 



Algebraic domains and double points 

A special case occurs when the Schwarz function on the physical sheet is 
meromorphic. It has no other singularities than poles of A. This is the case 
of algebraic domains [29]. They appear in the semiclassical case (see Sec. 3). 
This situation occurs if cuts on the physical sheet, situated outside physical 
droplets, shrink to points, i.e., two or more branching points merge. Then the 
physical sheet meets other sheets along cuts situated inside physical droplets 
only and also at some points on their exterior [double points) (see figure cap- 
tions for Figs. 3, 4). In this case the Riemann surface degenerates. The genus 
is given by the number of physical droplets only. The filling factors (40) are 
all positive. 

In the case of algebraic domains, the physical branch of the Schwarz function 
is a well-defined meromorphic function. Analytic continuations of z from dif- 
ferent disconnected parts of the boundary give the same result. In this case, 
the Schwarz function can be written through the Cauchy transform of the 
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Fig. 12. Critically degenerate Joukowsky map. 
physical droplets: 

S{z)^A{z) + l[J^. (46) 

TT Jd Z — Q 

Although algebraic domains occur in physical problems such as Laplacian 
growth, their semiclassical evolution is limited. Almost all algebraic domains 
will be broken in a growth process. Within a finite time (the area of the 
domain) they degenerate further into critical curves (Fig. 7, 12). The Gaussian 
potential (the Ginibre-Girko ensemble), which leads to a single droplet of the 
form of an ellipse is a known exception. 



4- 1-2 Critical degenerate curves 

Algebraic domains appear as a result of merging of simple branching points on 
the physical sheet. The double points are located outside physical droplets. Re- 
maining branching points belong to the interior of physical droplets. Initially, 
they survive in the degeneration process. However, as known in the theory 
of Laplacian growth, the process necessarily leads to a further degeneration. 
Sooner or later, at least one of the interior branching points merges with one of 
the double points in the exterior. Curves degenerated in this manner are called 
critical. For the genus one and three this degeneration is discussed below. The 
critical degenerations are depicted, respectively, in Fig. 12 and Fig. 7. 

Since interior branching points can only merge with exterior branching points 
on the boundary of the droplet, the boundary develops a cusp, characterized 
by a pair p, q of mutually prime integers. In local coordinates around such a 
cusp, the curve looks like ~ y^. 

The fact that the growth of algebraic domains always leads to critical curves 
is known in the theory of Laplacian growth (see e.g., [31]) as finite time sin- 
gularities. It is also known in the Hermitian one- and two-matrix models as 
critical points (intensively studied in 2-D-gravity [7]). 

The degeneration process seems to be a feature of the semiclassical approxi- 
mation. Curves treated beyond this approximation never degenerate. 
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4- 1-3 Simply-connected domains. Conformal maps. 

The case when the complement to the support of eigenvalues is a simply- 
connected algebraic domain in the extended (i.e. including oo) complex plane 
is particularly important. All the filling factors are zero except one, which is 
equal to A^, so there is only one droplet D. In this case, the Schottky double 
of the exterior of the droplet is a Riemann sphere. 

The exterior domain can be conformally and univalently mapped onto the 
exterior of the unit disc. We call this map w{z), and its inverse z{w). For 
algebraic domains, the inverse map is a rational function of w. Choosing the 
normalization such that z{oo) = oo, we represent the inverse map by a half- 
infinite Laurent polynomial 

z{w) — rw + UkW"'' , \w\ > 1. (47) 

fc>0 

The coefficient r, called the external conformal radius of the domain D, is 
chosen to be real positive. The Schwarz refiection of the inverse map is 

z{w-^) = rw-^ + \w\ > 1. (48) 

fe>0 

Given a point z, there are d values of w solving the equation (47), where 
d is the number of poles of the rational function z{w). Among them, only 
one solution corresponds to the conformal map. It is the solution such that 

w ^ z/r as z ^ oo. The global coordinate w on the Riemann sphere provides 
a uniformization of the Schottky double. In this coordinate, the involution 
reads w 1/w. The Schwarz function is 

S{z) = z{w~\z)). (49) 



Choose a regular point of the potential inside the droplet and expand 

V{z)^Ytk{z-^o)'. (50) 

fe>0 

The coefficients tk have a simple geometric interpretation. It follows from (44) 
that they are harmonic moments of the exterior of the droplet (with respect 
to the point ^o)- 




We also mention the area formula 

t^r^-Y.klu^'^l^- (52) 
fe>i 
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If the vector potential A{z) is a polynomial of degree d — 1, then the inverse 
map z{w) is a (Laurent) polynomial: u^'^'' = 0, k > d—l.li A{z) is a rational 
function with d — 1 simple poles in the finite part of the complex plane, the 
map z{w) is a rational function with d — 1 simple poles and one extra simple 
pole at infinity. 

4-2 The generating differential 
The meromorphic differential 

dn = S{z)dz (53) 

plays an important role. It is called generating differential [3,32]. According to 
(45), on the physical sheet it has the same poles and residues as the differential 
Adz: 

dn^ Adz + S4z)dz. (54) 

Below we use the following properties of the generating differential. 

(i) The periods over a-cycles (boundaries of the droplets) are purely imaginary, 
and are integer multiples of 27ri. They compute areas of the droplets (the 
filhng factors (40)): 

Zmn Jsia 

The filhng factors of physical droplets (belonging to the physical sheet) are 
positive. 

(ii) The real part of the integral of the differential {z — S{z))dz from some fixed 
point to a point on the boundary of a droplet (a point on a a-cycle) has 
the same value for all points of the boundary: 

4>a — + 2"^^ [ dVt — const, for all z e aQ,. 

This quantity does not depend on z, but does depend on unless is on 
the boundary. However, the difference 0q, — (pp depends on the a-cycles only. 
It is equal to a b-period of the differential dO,: 

(pa-(pp^ f dfl, (55) 

where ha/s is a cycle connecting aa and a^j cycles. 
For proofs and more details, see Ref. [3]. 

Periods over b-cycles (pa — <Po play a role of chemical potentials for the fill- 
ing factors. Here denotes a chosen reference droplet. One can use chemical 
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potentials to characterize evolution of the curve instead of filling factors. 



5 Examples 



5.1 Genus zero - ellipse 



The potential is Gaussian, V{z) = ^2-2^, 2|t2| < 1, A{z) = 2t2Z. It has a 
simple pole at infinity. The number of sheets of the Schwarz function is two. 
At infinity, -S"(^) is finite. Therefore, >S"(^) takes every value twice. It must 
have two zeros at most, i.e., there are two branching points. The Riemann- 
Hurwitz theorem says that genus is zero. The curve with two branching points 
and two sheets has the form z) = zz + kiz"^ + kiz"^ + /c2 = 0. Conditions 
that S{z) ~ 2t2Z + nhz~^ + . . . at 2; —> 00 for the physical branch and the 
unitarity condition (37) determine the curve in full: 



2 A-^\t' 
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fn{z, z) — zz — (toZ^ + toz'^) ; — uk [-^j- = 0. (56) 



The physical branch 



.<",)^....^^.(i-^i-p3^^) 

is the Schwarz function of ellipse. The second branch S^'^\z) (2t2)^^z does 
not correspond to any real curve. The droplet (see Fig. 2) is an ellipse with 
the quadrupole moment 2|t2| and area Trt — Trnh. 

Recurrence relations Eqs. (20, 18) truncate 

Zipn = Tnlpn+l + ^nV'n-l, (^'^V')n = rn-l^l'n-l + ^n+l^n+l- (57) 

Equations (19, 22) imply 

2^2 = ^^^^, nh ^ rl - lun+il"^ , (58) 



where the second equation is an analog of the area formula (52). Together 
they give a growth equation for the quantum analog of conformal radius nh — 
(1-|2^2|>^ 

The operator L"^ can be cast in the form of a 2 x 2 matrix (the number of 
sheets of the curve). Writing (57) for n+ 1, we express Xn-i and Xn+2 through 
X — iXn, Xn+i) from the first pair, and substituting them to the second pair 
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we obtain 



r-n. - 



Un+2Un+l 



Un+1 - 

Un+2 



Un 



(59) 



Formally, one can say that the recurrence relation (57) represents the shift 
operator as the 2x2 matrix 







Un Z 



The relation Cr 





With the help of (58), reads 



Un+1 



^ Un+2 



Wn then yields (59). 



z{2t2)-' -r„(2t2)-'(l-4|t2r 
r„(l-4|i2p) 2t2Z 



Computing the determinant det(£„(^) — z), we obtain the curve (56), already 
determined by the singularities. 

The recurrence relation can be used to generate the Hermite polynomials. 
The model with potential V{z) = t2Z^ + Qlog^; can also be solved explicitly 
(Ref. [39]). It generates the Laguerre polynomials. These are the cases where 
biorthogonal polynomials are classical orthogonal polynomials. 

Classical polynomials correspond to a genus zero curve. Higher genus curves 
generate more complicated, but more interesting polynomials. 

The classical limit of the recurrence relations gives a conformal map of the 
exterior of the unit disk to the exterior of the ellipse, z{w) — rw + — . The 



double- valued function Wi 2(2;) 



j_ 

2r 



Zi,2 = 



zi){z - Z2) 

becomes single- valued on a two-sheet covering of the 2;-sphere plane. The 
branch Wi{z) is such that u^i — > 00 as 2; — > 00. It defines the inverse map 
from the ellipse exterior onto the exterior of the unit disk in the w-plane. The 
function z{w~^) = rw~^ +uw is a mcromorphic function of w with one simple 
pole at -u; = 0. Treated as functions of the two sheets S^^''^\z) = z{wi^{z)) 
have two critical points 2:1,2 and constitute a sphere. Parameters of the map 
are related to the parameters of the potential through the classical hmit of 



(58): t2 



2r' 



t 



^2 _ 



u\^. The general formula for the semiclassical be- 



havior of the wave function (76) gives the known asymptotics of the Hermite 
polynomials [34]. 



The cut inside the droplet may shrink to a point. In this case only one (physi- 
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cal) sheet remains, and the Schwarz function becomes a rational function with 
one simple pole. The wave functions are monomials. The droplet is a disk. 

A more interesting degeneration occurs at the critical value of |t2| = |, when 
the two branching points reach the boundary of the droplet. The ellipse de- 
generates into the cut [zi, 2^2]. 

Beyond the critical value (|t2| > |), the ellipse appears again, but on the 
unphysical sheet. In this case, there is no physical droplet, but there is a 
virtual droplet. At this value of t2, the matrix integral diverges, but a solution 
of the differential equation (28) exists. It is not a polynomial, however. 

Below we use a more direct method to find the quantum curve bypassing the 
matrix representation of . The method reflects the triangular structure of 
, and closely resembles the procedure used in the classical case. We will use 
this method for the following, more complicated, examples. 

Apply (57) to an eigenvector c„ of the matrix Cn{z) with an eigenvalue z. The 
equations 

ZCn = Un+lCn-i-1 + T^n-lCn-l 

are compatible if the point {z, z) belongs to the curve. To obtain a compat- 
ibility condition, we express c„+i and c„_i in terms of c„. The results must 
differ by a shift n ^ n + 2 

This gives the curve (56). 



5.2 Genus one - an aircraft wing 

The potential is V{z) = — alog(l — z/ 13) — 72;, A{z) = — jz^ — 7- There is 
one pole &X, z — (5 on the first (physical) sheet. At 2; = 00 on the first sheet 
S{z) — > —7 -I- n^^. Therefore, the Schwarz function has another pole at the 
point —7 on another sheet. All the poles are simple. According to the general 
arguments of Sec. 4, the number of sheets is 2, the number of branching points 
is 4. The genus is 1. The curve has the form 

7(2;, z) — z^z'^ -\- kiz^z -\- kizz'^ + k2Z^ -\- ^2^^ + k^zz -\- k^z + k^z -\-h — Q. 



det 



Z Ur, 



det 



Tn+l Z 



— dndfi+li 



dn — det 
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The points at infinity and —7 belong to the second sheet of the algebraic 
covering. Summing up, 



S{z) 



a 


as 


z - 


-A, 




as 


z - 


OOl, 


nh—a 


as 


z - 


-72, 


(/9-f) 


as 


z - 


002. 



where, by 1 and 2 we indicate the sheets. 

Poles and residues of the Schwarz function determine all the coefficients of the 
curve f{z, z) = a{z){z - S^'-\z)){z - S^''\z)) = '^{z - - S^''\z)) 

except one. The behavior at 00 of 2;. ^ gives fci = 7 — /9, /c2 = —7/?. Hereafter 
we choose the origin by setting 7 = 0. The equation of the curve then reads 
fn{z, z) = 0, where fn{z, z) is given by 

z^z'^-z'^z^-zz'^(3+{\^^\'^ + a + a - n?i) zz+z^{nh-a)+z(3{nh-a)+hn (60) 

The free term is to be determined by filling factors of the two droplets Vi 
and U2 = n — ui. K detailed analysis shows that the droplets belong to different 
sheets (Fig. 10). Therefore, 1^2 is negative. 

A boundary of a physical droplet is given by the equation z = S^^\z) (Fig. 3). 
The second droplet belongs to the unphysical sheet. Its boundary is given by 
z = S^'^\z). The explicit form of both branches is 

c(i 2) _ 1 5 ^(""^ - a) + (a + a - nh)z ^ ^J{z- zi){z - Z2){z - Z3){z - Z4) 
~r 2{z-(3)z ' 

where the branching points Zi depend on 

If the filling factor of the physical droplet is equal to n, the cut inside the 
unphysical droplet is of the order of v^- Although it never vanishes, it shrinks 
to a double point z^ = = z^. in a semiclassical limit. The sheets meet at the 



double point rather than along the cut: y {z — Zi) {z — Z2) {z — z^) {z — Z4) 

(z — z^)^ [z — Zi){z — Z2). In this case, genus of the curve reduces to zero 
and the exterior of the physical droplet becomes an algebraic domain. This 
condition determines /i, and also the position of the double point (Fig. 11). 
The double point is a saddle point for the level curves of f{z,z). If all the 
parameters are real, the double point is stable in x-direction and unstable in 
y-direction. A saddle point is a signature of a virtual droplet [4]. 

If this solution is chosen, the exterior of the physical droplet can be mapped 
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to the exterior of the unit disk by the Joukowsky map 



z{w) = rw + Uq + 



u 



w — a 



\w\ > 1, \a\ < 1. 



(61) 



The inverse map is given by the branch wi{z) (such that — > oo as z — > oo) 
of the double valued function 



z - Uo + ar ± ^J{z - Zi){z - Z2) , ^1,2 = Uo+ar^2^r{u + auo) 



The function 



z{'w ^) — rw ^ + Uo + 



u 



(62) 



— a 

is a meromorphic function of w with two simple poles at w = and w — 
a~^. Treated as a function of z, it covers the 2;-plane twice. Two branches 
of the Schwarz function are S^^''^\z) = z{wil{z)). On the physical sheet, 
S^^\z) = z{wi{z)) is the analytic continuation of z away from the boundary. 
This function is meromorphic outside the droplet. Apart from a cut between 
the branching points Zi^2, the sheets also meet at the double point z^ = —7 + 
Q-i^g2i0^ where S^^\z^j = S^'^^z^), = arg(ar + jz^)- 

Analyzing singularities of the Schwarz function, one connects parameters of 
the conformal map with the deformation parameters: 



7 

nh — a 



2 _ HI 
r _|_ _|_ ua 



(63) 



Area of the droplet ^ nh — — 





u 


2 


(1- 




a 


2^2 



A critical degeneration occurs when the double point merges with a branching 
point located inside the droplet (2;* = Z2) to form a triple point 2;^,^. This may 
happen on the boundary only. At this point, the boundary has a (2, 3) cusp 
(Fig. 12). In local coordinates, it is ~ y^. This is a critical point of the 
conformal map: w'{z^^) — 00. A critical point inevitably results from the 
evolution at some finite critical area. 

A direct way to obtain the complex curve from the conformal map is the 
following. First, rewrite (61) and (62) as 

z — Uq -\- ar — rw -\- a{z + 7)^/;"^ 
< (64) 

z — Uq + ar — rw~^ + a{z + 7)^, 
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and treat w and 1/w as independent variables. Then impose the condition 



w ■ w 



1. One obtains 



det 



z — uq + ar a{z + 'j) 
z — UQ + ar r 





2 


/ 






det 









r a{z + 7) 
a{z + 7) r 



This gives the equation of the curve and in particular /i, in terms of Uq, r, a 
and eventually through the deformation parameters a, /3, 7 and t. 

The semiclassical analysis gives a guidance for the form of the recurrence rela- 
tions. Let us use an ansatz for the L-operator, which resembles the conformal 
map (61): 

L = rnW + + {W - an)~^Un, 

so that 



(w - a„)L = {w - an)rnW + (w - an)u^n^ + Un, 
where w is the shift operator n — > n + 1. 



(65) 
(66) 



Now we follow the procedure of the previous section. Since the potential has 
only one pole, can be cast into 2x2 matrix form. Let us apply the lines 
(65, 66) to an eigenvector (c„,c„_|_i) of a yet unknown operator jC„, and set 
the eigenvalue to be z: 



{z + r„_ia„_i - 'ui°))c„ = r„c„+i + a„_i(^ + 7„_i)c„_i 



(z + r^an - xii+i)c„ = a„+i(z + 7„+i)c„+i + r„c„_i. 



(67) 



We have defined 7n = ^ — u^- The equations are compatible if c„_i and c„+i 
found through c„ differ by the shift n — > n + 2. We have 



-71+ 1 



;idet 



z + r„_ian_i - a„_i(2; + 7„_i 
~ , - _ -(0) 



Cn-1 = det 

On 



r„ 2; + Tn-ian-l - 

a„+i (5 + 7„+i) z + r„a„ - M^fji 



(0) 



On 



--n J ) 



where 



dn — det 
This yields the curve 



r„ a„_i(2; + 7„_i) 

a„+i(^ + 7„+i) r„. 



(68) 
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Comparing the two forms of the curve (60) and (68), we obtain the conserva- 
tion laws of growth: 



7 = 7n = 3 u. 



ni 



(0) 

o _ I (0) '^n+l'^n^n+l 

P — ~ r l^n+l + ~. = ) 

nti-a^ r„rn+i . 

They are the quantum version of (63). 

The poles of the vector potential field determine most elements of the matrix 
Cn- We have 



where A, B and C are ^-independent 2x2 matrices. Their elements are An = 
p, A21 = pan, A12 = A22 = 0, Bu=t- a, Bi2 = B21 = B22 = 

_ tin — 1 

0, trC = —a, detC = 0, Tin = —Cn — j3an{Bi2 + C12). The undetermined 
matrix elements are found similarly to the case of the ellipse. We treat the 
recurrence relation (the first fine of (67)) as a 2 x 2 representation of the shift 
operator 

\ 



1 

^ an 



and compute Cn from the second line of (67). The matrix elements of C„ 
appear to be rather cumbersome. 



5.3 Genus three - hypotrochoid 



The potential is V(z) = t^z^. A general polynomial potential, including this 
example, has been analyzed in Ref. [9]. We briefly review it here to illustrate 
features of the curve with multiple poles of the Schwarz function. In this case, 
A{z) has a double pole at infinity. The number of sheets is three. The behavior 
of the Schwarz function at infinity is S{z) ~ St^z^ and the symmetry of the 
potential under 27r/3-rotation restricts the equation of the curve to the form 

fn{z, Z) = {zzf ~ ^ ~ ^ + ^nZZ + K ^ 0, (69) 

with two unknown coefficients. Since S{z) has a double pole at infinity on the 
physical sheet, it has a branching point at infinity on the other sheets. This is 
the point where two unphysical sheets meet. Eq.(69) suggests that there are 
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9 more branching points. Therefore, genus of the curve is g = 3. There are 
four droplets, but only one belongs to the physical sheet. 27r/3 rotations leave 
the physical droplet invariant and exchange unphysical droplets. Therefore, 
the filling factors of unphysical droplets must be equal. This number, together 
with the filling factor of the physical droplet, determine the coefficients in 
(69). 

Let us start from the semiclassical analysis. Assume that the filling factor of 
the physical droplet is n. We then have an algebraic domain bounded by a 
hypotrochoid. Fig. 4. Six critical points sitting outside the droplet collapse 
two by two to three double points. The function 

u 

z{w) = rw -\ (70) 

maps the exterior of the unit disk to the exterior of the hypotrochoid. The 
inverse function has three branches Wi{z), i = 1,2,3. At infinity, the three 
sheets have the leading terms Wi{z) z/r, W2^3{z) i^u/zY^"^. Therefore, 
the inverse of the map (70) is Wi{z). The branch S{z) = z{wi^{z)) is the 
Schwarz function. It is a meromorphic function outside the domain with one 
double pole at infinity, S{z) ~ 3^32;^ + {nh)/z. This gives u — Stsr"^, while the 
area is = — 2|iip = — ISj^apr^. The coefficients of the equation of the 
curve are determined by 

, (l-9|t3|V)(l + 18|t3|V) r^(l-9|t3|V)3 

~ nli 19 ' nli 19 ■ \ ) 



Four remaining branching points are z = 00 and three points solving dz = 
(r — 2uw~^)dw = 0, at = 2u/r. They are situated inside the droplet. Three 
double points solve S^'^\z) — S^\z). They are outside the droplet: 



(r2 




\u 


2^3 




u\ 





At a critical degeneration the double points merge on the boundary to the 
branching point inside the droplet, z^ — z{w^), = 1 (Fig. 7). 

Quantum analysis starts from recurrence relations guided by the classical case. 
Eqs. (20, 18) read 

Zll^n = rni>n+l + Uni>n-2: {L^i>)n = Tn-ll/jn-l + Un+2i>n+2- (72) 

Acting on the eigenvector of we get 
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Writing the first equation for n — > n ± 1, we express c„+2 and Cn-i through c„ 
and Cn+i to obtain 

ZCn — Z Cn+l i [Tn-lTn+l — 'U'n+lUn+2 )Cn-l, 

fn+1 fn+1 



V rn-2 J rn-2 

Equations are compatible if z belongs to the curve 



det 



( r ^l^+l^w+2 ^ 



where 



r„+2 



dn — det 



■ det 



rn-2 



V rn+1 



V " rn-2 ) 



rn-2 



Comparing with (69), we find the quantum version of (71): 



ZUn+2 
rn+1 



l-9|i3r(<-i + <+i) + 



1 - 9|t3iX-i)(i - 9|t3r<)(i - gitsiv 



12^2 N 
n+lJi 



where we used the conservation law and a quantum analog of the area formula 

Un+l 



3U 



hn. 



Together, these equations give a closed equation for r^, 
which is a discrete analog of the Painleve I equation. 

The matrix form of can be obtained in a similar manner. It is a 3 x 3 matrix 
acting on the 3- vector {xn,Xn+i,Xn+2y- Its analytic structure is 

and the matrices A, B, C are 
/ 









r„_i(l - \3h\^rl) 

r„(l - \3ts\'rl^,) 








35 



^(3^"3) 

|3i3|V„r„_i 








'^0 ^ 




3t3 



6 Semicleissical wave function 



In the scmiclassical limit, h ^ 0, while t = hN, t^°'^ = hua, and the potential 
V{z) are kept fixed. The quantum curve and evolution equations go to their 
semiclassical counterparts. Also, the semiclassical wave functions are written 
through the objects of the semiclassical theory, i.e., through differentials on 
the semiclassical curve. 

A technical difficulty is that wave functions with close indices, say tpn and 
ipn+i, are not necessarily close to each other as n — > oo. Choosing a sequence 
of states which has a semiclassical limit is sensitive to the potential, and refiects 
the configuration of semiclassical droplets. This problem does not arise in the 
case of a single droplet in the algebraic case, where states with close numbers 
are close. As we discussed above, this case corresponds to degenerate Riemann 
surfaces. Hereafter, we assume that states with close n's are close, and proceed 
to the semiclassical limit. The result obtained for this case can be generalized 
to any smooth Riemann surface. 

First, we rewrite the L operator (20) in a more symmetric form L — w~^^'^rw 

1/2+ 

Y^^-k/2^{k)^-k/2 ^ by redefining r and u^^\ Then we search for a semiclassical 
expansion 

(in what follows, we suppress the index n). 

In the first two leading orders, the L-operator (20) reads 

Ti 

L — > z{w) + -{wdy,z{w), dt]Ai + . . . , 

where logtu = dtAo and { , } means anticommutator. The function 

z{w) = rw + ^ u'^'>w~^ (73) 

fc>0 

is obtained by replacing the shift operator by its classical value w — > w. 
The recurrence relation Lip{z) — zip{z) in two leading orders reads 

z{w) = z, dtAi = --^ — logw = dtAo. 

2 a„ziw) 
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It gives 

A,{z) = \\og w'{z), (74) 
where w{z) is a multivalued function inverse to (73). 

In the leading orders, (18) reads 

hd,xl^{z) = ( - I + Kw-^) + \{wd^z{w-^). dt} + A, + ... )V'(^). (75) 

The function z{w'^) is the Schwarz reflection of z{w) with respect to the 
contour whose exterior is mapped to the exterior of the unit disk by the func- 
tion 'w{z). This map is discussed in Sec. 4.1.3. The contour is the boundary 
of a semiclassical droplet. We observe that in the simplest semiclassical limit 
the conformal map and its Schwarz reflection arc limits of the L and L"!" op- 
erators of Sec. 2. The complex curve in the semiclassical limit is given by 
f{z,z) = a{z)Y[{z — z{w~^{z)), where the product is taken over all possible 
branches of the multivalued function w{z). 

Introducing the Schwarz function S{z) = z{w~^{z)), as in (49), we obtain 

Ao = -^ + n{z), 

where Q,{z) — J^^ dQ is the primitive function of the generating differential 
introduced in Sec. 4.2. The point is such that log w(^o) — 0. It is a boundary 
point. The integral goes along the physical sheet. 

Summing up, we obtain the result: 

^|J{z) ~ (27r3;i)-Vy^^'(^)eH(-il^l'+^W). (76) 

Formula (76) was reported in [2]. A similar formula for the two- matrix model 

was obtained in [33]. Being applied to the Hermite or Laguerre polynomials 
(sec example in Sec. 5.1), the formula gives the asymptotical behaviour found 
by Tricomi in 1941 [34]. 

The overall numerical factor can be fixed cither by comparison with asymp- 
totics of biorthogonal polynomials for the Ginibrc ensemble, or through nor- 
malization of the semiclassical wave function. The latter is described below. 

An extension of this formula to the case of a general complex curve and multi- 
ple droplets is straightforward. The function Q{z) should be again understood 
as an integral of the generating differential dQ. The integration path starts at 
a point on the boundary of some "reference" droplet and ends at the point z. 
The correction factor, w'{z) in (76), is replaced by the meromorphic Abelian 
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differential of the third kind W^°°'°°\z), introduced in Sec. 4.2. As a result, 
we obtain: 

^|J{z)^/d^ ~ ^iy(oo,db) (^) g-i i^-Iio ^^'^'^') . (77) 

By virtue of (40), monodromies of this wave function around droplets are 
g27rii^a , where are fiUing factors. 

The asymptotical representation (77) fails around cuts or double points. It 
is valid only around boundaries of physical droplets and only if the distance 
between boundary of a physical droplet and a cut or a double point is much 
larger than ^/h. In other words, it is valid for smooth non-critical curves. The 
scmiclassical asymptotics of the wave functon in the whole plane is far beyond 
the scope of this paper. 

The following properties of visualize the result. 

(i) The real part of is constant along boundary of any droplet (see Sec. 4.2), 
generally different for different droplets 

2neAo{z) ^(j>a, ze dD^. 

(ii) Derivatives and of TZe Aq vanish on the boundary: 

d,TZeAo{z) = d^TZeAo{z) = 0, ze dD. 

This is the condition for scmiclassical extrema (41). 

(iii) In the leading order, the change of TZe Aq under small deviations along the 
normal direction away from the boundary to both exterior and interior does 
not depend on the point of the boundary: 

d'^TZeAoiz) = -2. 

(iv) The higher order expansion away from the boundary depends of the curva- 
ture of the boundary k{z),z e dD^- 

TZeAoiz + 6n) = 0c. - {Snf + h{z){Snf - ^k\z){dn)^ + . . . , (78) 

where 5n is a normal deviation from the boundary directed outward. This 
expansion fails at very curved parts of the contour, where k{z) ~ h~^^'^. 

Altogether, these properties mean that the wave function peaks at the bound- 
ary of a physical droplet and decays as Gaussian at the rate \/h (see Fig. 6): 

\iP{z+5n)\'' r^e^e-i^^-^\ z e dD^. (79) 
The amplitude of the wave function may be different on different boundaries. 
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Since the wave function is localized on the boundary, one can integrate in the 
normal direction. In the case of one droplet, we have 

J mz)\'d{6n) = ^dnlog\w{z)\ = ^\w'{z)l z e dD (80) 

(here 9„ is the normal derivative). Further integration along the boundary 
is consistent with the normalization / \t/jN{z)\'^(Pz — 1. This can be used to 
determine the constant in (76). 

According to the arguments of Sec. 1.6, Eq. (80) can be interpreted as a rate 
of growth. 

Note that Eq. (74) and the next to leading order of (75) give the relation 

w{dyjzdtz - dtzdyjz) = 2, (81) 

where one differentiates at a fixed w. This relation represents evolution of the 
growing domain as a time-dependent conformal map. The area of the domain 
grows linearly with time, while its harmonic moments stay fixed. This type of 
classical growth is equivalent to the Laplacian growth. Relation (81) has been 
known in the literature for a long time [38]. It is a classical hmit of Eq. (21). 



7 Appendix. Laplacian growth 



Laplacian growth referres to growth of a planar domain, whose boundary prop- 
agates with a velocity proportional to the gradient of a harmonic field. The 
Hele-Shaw problem is a typical example (see [31] for a review). The exper- 
imental set-up consists (Fig. 13) of two parallel and horizontal glass plates, 
separated by a small distance along the vertical direction, called the gap, b. 
The gap is small enough compared to typical intermolecular distances in fluids, 
so that the system can be considered two-dimensional. 

Initially, the space between the plates is filled with fluid 2, of large viscosity 
(tipical examples are silicon oil or liquid crystals). Through a small opening 
in the middle of the upper glass plate, fluid 1, of much lower viscosity (such 
as air), is inserted into the system, at a constant rate. Both fluids are incom- 
pressible and immiscible, so fluid 2 must be evacuated from the system at the 
same rate as fluid 1 is being introduced. This is done somewhere at the edge 
of the glass plates. 

Given that the fluids are incompressible and immiscible, the dynamics of the 
system is reduced to that of the interface between them. Depending on the 
relevant parameters of the experiment (the coefficients of viscosity, /Xi^2 of the 
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Fig. 13. The Hele-Shaw setup. 




Fig. 14. Laminar Hele-Shaw dynamics. 

fluids, the coefficient of surface tension cr, the gap b and the pumping rate Q), 
the interface dynamics wih exhibit very different behaviors. At low pumping 
rates Q, for given cr, the interface will remain smooth at all times, regardless of 
how much fluid is being inserted (Fig. 14). However, at high Q (or equivalently, 
at very low a), the boundary will develop the so-called fingering instability, 
where certain regions grow in the shape of fingers, whose tips split and form 
secondary fingers, etc. Real experiments (Fig. 5) show that in this regime, the 
interface seems to evolve into a self-similar, fractal structure. 

Under the assumption that the system is two dimensional, the equations gov- 
erning the dynamics of the interface lead to Darcy's Law: the velocity field Vi 
and pressure Pi in fiuid i are related as: 

b^ ^ 

Vp,, (82) 



and therefore the continuity equation accross the interface leads to 

v-n^-——n-Vpi^-——n-Vp2, (83) 

where n is the unit vector normal to the interface at a given point. 

In each fiuid, the pressure field is supposed to solve the Laplace equation, 

Api = 0, (84) 
with appropriate boundary conditions. For instance, since fiuid 2 is drained 
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from the system somewhere away from the origin, fluid 2 has a sink at infinity, 
which leads to the condition 



^2(2;) ^ -log |2;|, z ^ 00. 



(85) 



The coefficient of surface tension enters the problem through the Laplace for- 
mula relating the pressure jump accross the boundary and the local curvature 
k{z): 



The first simplification of the problem consists in setting the viscosity of fluid 
1 to zero, fii = 0. The continuity condition then implies pi — constant, and 
we may therefore choose pi = (a redefinition of ^2)- 

The second important simplification is neglecting the effect of surface tension. 
This assumption has important physical meaning, since the resulting prob- 
lem should necessarily lead to turbulent dynamics. Mathematically, the two 
assumptions lead to the following reformulation of the problem: 



Ap = on D_, p = 0, K ~ dnP on F, p{z) — log \z\, (87) 



where Vn is the normal component of the velocity, is the domain occupied 
by fluid 2, and F is the boundary between fluids. Under these assumptions, in 
order to describe the interface dynamics, we must solve an exterior Dirichlet 
problem for the pressure field. 



Acknowledgments 

We are indebted to A. Kapaev, V. Kazakov, I. Krichever, I. Kostov, A. Mar- 
shakov and M. Mineev-Weinstein for useful discussions, interest in the subject 
and help. P.W. and R.T. were supported by the NSF MRSEC Program under 
DMR-0213745, NSF DMR-0220198 and by the Humboldt foundation. A.Z. 
and P.W. acknowledge support by the LDRD project 20020006ER "Unstable 
Fluid/Fluid Interfaces" at Los Alamos National Laboratory and M. Mineev- 
Weinstein for the hospitality in Los Alamos. A.Z. was also supported in part 
by RFBR grant 03-02-17373 and by the grant for support of scientific schools 
NSh-1999.2003.2. P.W. is grateful to K.B. Efetov for the hospitality in Ruhr- 
Universitaet Bochum and to A.Cappelli for the hospitality in the University of 
Florence, where this work was completed. We are grateful to Harry Swinney 
for permitting us to use Fig. 5 from [16]. 



P2{z) -Pi{z) = aK{z). 



(86) 
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